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Abstract 

The minimal energy B = 6 solution of the Skyrme model is a static soliton with D4,d 
symmetry. The symmetries of the solution imply that the quantum numbers of the ground 
state are the same as those of the Lithium-6 nucleus. This identification is considered further by 
obtaining expressions for the mean charge radius and quadrupole moment, dependent only on 
the Skyrme model parameters e (a dimensionless constant) and _FV (the pion decay constant). 
The optimal values of these parameters have often been deliberated upon, and we propose, 
for B > 2, changing them from those which are most commonly accepted. We obtain specific 
values for these parameters for B = 6, by matching with properties of the Lithium-6 nucleus. 
We find further support for the new values by reconsidering the Q-particle and deuteron as 
quantized _B = 4 and B = 2 Skyrmions. 



1 Introduction 

A nonlinear theory was put forward by T. H. R. Skyrme in 1961 [21] • It sets out to describe a low- 
energy effective theory of QCD, and the Lagrangian of the theory is defined in terms of pion fields. 
The most successful and promising physical interpretation of the theory is as a description of atomic 
nuclei. The Skyrme model admits topological soliton solutions, falling into sectors labelled by an 
integer-valued topological degree. This quantity is identified with baryon number. A quantized 
Skyrmion of topological charge B is interpreted as a nucleus with baryon number B. 

To date, relatively little work has been done on the electromagnetic properties of quantized 
Skyrmions and their comparison to known experimental results. Two papers by Braaten and 
Carson [HJ Ej consider the electrostatic properties of the B = 2 Skyrmion, and their comparison 
to the experimentally determined properties of the deuteron. This analysis was extended in |19j 
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to account for the softest vibrational modes of the B — 2 Skyrmion, the quantization of which 
resulted in an accurate prediction of the mean charge radius of the deuteron. 

In this paper we perform a similar analysis for the B — 6 Skyrmion. However, performing a 
quantization of all possible vibrational modes for _B > 2 is technically very difficult. We therefore 
consider rescaling the Skyrme model parameters (introduced in the next section), whilst sup- 
pressing the vibrational modes, in order to calculate observables that more closely correspond to 
experiment. 

The outline of the paper is as follows. In the following section we present the Skyrme model and, 
following give an expression for the electromagnetic current. Section 13 provides an overview 
of the rational map ansatz for Skyrmions ^51, which is used in later sections. Section ^discusses 
the fermionic quantization of solitons, as applied to the Skyrme model, recalls the Finkelstein- 
Rubinstein (FR) constraints jj^l i and reviews Krusch's method of determining the PR constraints 
|17| . In sectional we recall the semi-classical method of soliton quantization, following (2), in which 
only the translational, rotational, and isorotational collective coordinates are considered. In section 
El this quantization of the B = 6 Skyrmion is performed, and we find that the symmetries of the 
minimal-energy classical solution, via the FR constraints, ensure that the lowest-lying quantized 
state has spin 1 and isospin 0, the quantum numbers of the Lithium-6 nucleus. In section [7| we 
present an expression for the mean charge radius of the quantized B — 6 Skyrmion. Section 
IHl presents an expression for the quadrupole moment of the quantized B = 6 Skyrmion. The 
reparametrisation of the Skyrme model in the B = 6 sector is considered in sectionlHl and the mean 
charge radius and quadrupole moment of the quantized B — 6 Skyrmion are explicitly calculated 
and compared to experiment. In section [TUI we reconsider the B = 2 and _B = 4 Skyrmions in the 
light of the reparametrisation. We provide a conclusion in section ^2 The calculations presented 
in this paper support our interpretation of the Lithium-6 nucleus as the lowest-lying quantum state 
of the B = 6 Skyrmion. They confirm that if one wishes to apply the Skyrme model to nuclei, 
then a reparametrisation is desirable. 



2 The Skyrme Model 

The Skyrme model is a nonlinear theory of pions in three spatial dimensions, defined in terms 
of an S'J7(2)-valued scalar, the Skyrme field 1251 122| . It provides a low energy effective theory of 
QCD, becoming exact as the number of quark colours becomes large [271 128) . Its topological soliton 
solutions, Skyrmions, can be interpreted as nucleons and nuclei. 

The SU{2) Skyrme model has the Lagrangian density 

£ = § Tr a^C/9^[/-i-H^ Tr [d^UU-\d,UU~%d^UU-\ d''UU-^] + ^mlF-^ Tr ([/-I2) , (1) 

where U{t,x) is the Skyrme field, F^^ is the pion decay constant, e is a dimensionless parameter 
and is the pion mass. 

It is helpful to use energy and length units of F^/4e and 2/ei^7r respectively. In terms of these 
Skyrme units the Skyrme Lagrangian becomes 

J l^-^TT{R^Rn + ^TT{[R^,R,][R^,R''])+m^Tr{U~l2)^d^x, (2) 

where we have introduced the 5u(2)-valued current i?^ = {d^U)U^^ , and defined m = 2m7r/eF^. 

The scalar field U, at a fixed time, is a map from R'^ into S^, the group manifold of SU{2). 
However, the boundary condition t/ ^ I2 at spatial infinity implies a one-point compactification 
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of space, so that topologically U : ^ S^, where the domain is identified with U {00} . 
Thus, configurations U obeying this boundary condition fall into topological sectors labelled by 
their topological degree 

B = J BQ{x)d^x (3) 

where 

B,,{x) = ^ e^^^iiTTd-'UU-'d^UU-'d^UU-' . (4) 

The degree B, which takes integer values, is identified with baryon number. We refer to Bq as the 
baryon density. The minimal energy, static solutions of the field equation for each nonzero baryon 
number we call Skyrmions. 

The internal symmetry of the Skyrme model is the global isospin symmetry U — > AUA^ where 
A is an SU{2) matrix. This is generated by the infinitesimal transformations U U + ie.[Tp, U], 
where Tp : p = 1, 2, 3 are the Pauli matrices. To couple electromagnetism to the Skyrme model, the 
U{1) symmetry generated by the third component of isospin is gauged, with the derivatives in |^ 
and (0) replaced by the covariant derivative I?^f/ = df^U — ieoAfj,[Q, U], where Q = 5T3 + is 
the charge matrix of quarks. In the presence of electromagnetism, the form of the baryon current 
density given previously is unsatisfactory as it is not gauge invariant under 

U{x)~^U{x)+ieoa{x)[Q,U], -> + 9,,Q! . (5) 

The gauge invariant, conserved generalisation is given by |H] 

BJx) ^ -^ef,^afilT^i-U-^d''UU-^d°'UU-^d'^U + 3ieoA''TTQ{U'^d°'UU-^d'^U 

24:71 L 

-d°'UU-^d'^UU-^) + Sieod'^A" TiQiU^^d'^U + a^t/C/^^)} . (6) 
The electromagnetic current is determined as 

+ ^£^.0/35"^^" Tr {Q2a%[/-i + Q^U-^d^U + ^^Qd^UQU-^ - ^^QUQd^U-^] + ll{x) . (7) 



We are interested in static or slowly varying Skyrme fields, quantized to have isospin zero. In 
this case the electromagnetic current simplifies. It can be shown, assuming Iq(x) = 0, Aq ^ 0, 
A = and OqA^ = 0, that 

" " ^eo.yfeep,,7f tJtI: = ^ det r , (8) 
where the matrix V has elements 7f defined by 

(9,C/)C/-i = -i7frp . (9) 
Therefore, the electric charge density is half the original baryon density. 



3 The Rational Map Ansatz 

The rational map ansatz |14) applies the topological notion of suspension, using a rational map 
from S"^ S'^ to construct an approximate Skyrmion, a map from M.^ — s- S^. It exhibits a 
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nonlinear separation of variables, separating the angular and radial dependence of the Skyrme 
field. One identifies the domain S'^ with concentric spheres in R^, and the target S'^ with spheres 
of latitude on S''^. A point in R'^ can be denoted by its coordinates (r, z), where r is the radial 
distance from the origin and z specifies the direction from the origin, a point on the unit sphere. 
Via stereographic projection, the complex coordinate z can be identified with conventional polar 
coordinates by z = tan(0/2)e"^. Equivalently, the point z corresponds to the unit vector 



1 + \Z[' 



;{Z + Z,i{z- Z),l- \z\^), 



and conversely 



(n^)i + i{n,_)2 
1 + (n.)3 



(10) 



(11) 



The ansatz for the Skyrme field depends on a rational map R{z) — p{z)/q(z), where p and q are 
polynomials in z, and a radial profile function /(r). The value of the rational map R is associated 
with the unit vector 

nR^^-^^iR + R,i{R-R),l-\R\'). (12) 

The ansatz is 

U{r,z] 

with f{r) satisfying /(O) = vr and /(oo) 



: exp (i/(r)nij(^) • r) , 
0. 



(13) 



The degree of the rational map, N, is the greater of the algebraic degrees of the polynomials 
p and q. N is also equal to the topological degree of the map R viewed as a map from 5^ ^ 5*^. 
The baryon density of the ansatz is given by 



-/' /sin / 



and so the baryon number is given by 



B 



27r2 



1 + 



dR 



dR 



dz 



dz 



2i dz dz 



r dr . 



where 2i dz dz / (1 + \z\'^)'^ is equivalent to the usual 2-sphere area element sin 9 d9 1 
forward to show that B — N . 



(14) 
(15) 

It is straight- 



An SU{2) Mobius transformation Mi of R corresponds to an isorotation; an SU{2) Mobius 
transformation M2 of z corresponds to a rotation in physical space. They induce the following 
transformation of the rational map i?(z): 



R{z) ^ R{z) = AhiRiA'hiz))) . 



(16) 



The energy for a field of the form H13|l is 



E : 



which can be simplified to 



\R\' 



dR 



dz 



sin* / / 1 + |z| 



+2to2(1 - cos/) 



dR 



1 + 

2i dz dz 



dz 



dr , 



E = 4tt I [r^f^ + 2B sin^ /(/'^ + 1) + X + 2mV^ (1 - cos f)]dr. 



(17) 



(18) 
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Here T denotes the purely angular integral 



1 

47r 



1 + 



dR 



dz 



2i dz dz 

(l + |z|2)2 



(19) 



which only depends on the rational map R. Both B and X, and hence the energy E, are invariant 
under the transformations 

To minimise E one first minimises 2 over all maps of degree B. The profile function / may then 
be found by solving the second order ODE that is the Euler-Lagrange equation for the expression 
(|18|l with B and T as fixed parameters: 



(r2 + 2B sin^ /)/" + 2rf + sin 2/ B{f'^ - I) - I 



sin^ / 



— rri^r^ sin / = . 



(20) 



Note that the inclusion of the pion mass term in the Lagrangian density has no effect on the 
rational map, but results in the profile function being slightly modified, leading to higher energies 
than in the massless case. The optimised rational map ansatz has been shown to give a good 
approximation to the true Skyrmion for baryon numbers up to i? = 7 (and for a larger range of B 
in the massless pion case) 4, 14 . For B = 8,9 and beyond, and m of order 1, the structure of the 
minimal energy Skyrmions differs qualitatively from that given by the rational map ansatz pill3|. 



4 Fermionic Quantization of the Skyrme Model 

We recall that the configuration space C of the Skyrme model has connected components Cb, which 
are the homotopy classes, labelled by baryon number. The fundamental group of each component 
satisfies ttiICb) = ^2, and therefore all Cb admit double- valued functions. However, these double- 
valued functions can be defined as single- valued functions on Cb, the universal covering space of 
Cb- 

We write q, q' for different points oi Cb covering the point q & Cb- The condition that q ^ q' 
implies that a path from q to q' projects to a non-contractible loop in C^. Double- valued functions 
on Cb are thought of as functions on C^: 

^■.Cb^C, *-*(g). (21) 

Fermionic quantization requires that the wavefunction is defined on Cs and satisfies 

^(q) ^ ~^iq') . (22) 

Let a rational map R have the following symmetry, for some particular Afi and M2, 

R{z) ^ Kh{R{M2{z))) , (23) 

where Mi corresponds to an isorotation by 9i around rii and M2 corresponds to a rotation by 62 
around n2. For 62 ^ ink for k £ Z, M2 only leaves the antipodal points 

1 + (n2)3 1 - (n2)3 

fixed. Similarly, Mi only leaves i?±ni fixed, where -R±ni are defined similarly as above. Therefore, 
for the symmetry R{z) — Mi{R{M2(z))) to hold, we have 

i?(z_„J = R±n, - (25) 
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By reversing the signs of rii and 9i if necessary, we may set 

-^(^-na) = ^ni • 



(26) 



Now consider the Skyrme field configuration defined by the rational map ansatz, using this 
rational map and some profile function /, and assume the quantum wavefunction is non-vanishing 
for this field configuration. The symmetry gives rise to a loop in configuration space (one thinks of 
this as a loop in Cb by letting the isorotation angle increase from to 0i, with the rotation angle 
increasing from to 02) ■ This leads to the following constraint on the wavefunction: 

e^e2n2■■L^^e^n,■■K^^^ = XFi?|^') , (27) 

where L and K are the body-fixed spin and isospin operators respectively, and the Finkelstein- 
Rubinstein (FR) phase xfr enforces the fermionic quantization condition: 

_ ( +1 if the loop induced by the symmetry is contractible, , , 

= \ -1 otherwise. ^^^^ 

We note the following general result, proved in |17| : 

The value of xfr for a given symmetry of a rational map only depends on the isorotation angle 
9i and the rotation angle 02, where the angles are defined such that i?(z_n2) = -^ni, and is given 

by 

XFR ^{-if, where = ^{802 - 0i) . (29) 



5 Semi- Classical Collective Coordinate Quantization 

The Skyrme Lagrangian is invariant under the Poincare group of (3 -I- l)-dimensional space, 5'0(3) 
isorotations and some discrete parity transformations, which will not be considered here. Simi- 
larly, the space of static solutions, that is, configurations which minimise the energy functional, 
is invariant under the Euclidean group and isorotations, E3 x 50(3)^. By acting with the latter 
symmetry group on a static Skyrmion Uq we generate a set of new solutions: 

C/(x) = ^i(7o(I?(A2)(x- X))4 , (30) 

where Ai,A2 are SU{2) matrices and D{A2)ij — ^TT{TiA2TjA2) is the associated 5*0(3) rotation. 
The classical degeneracy of the above solutions is removed when the theory is quantized. We think 
of the parameters X(t), Ai(t) and A2{t) as dynamical variables and then quantize the resulting 
dynamical system according to standard canonical methods 0. This is the collective coordinate 
quantization of the Skyrmion Uq. As we shall only be concerned with the computation of spin 
and isospin, we shall ignore the translational degrees of freedom and quantize the solitons in their 
zero-momentum frame. We shall also ignore possible deformations of the Skyrmion Uq, which lead 
to vibrational excitations. 

Our dynamical ansatz is then U{x,t) = Ai{t)Uo{D{A2{t))x)Ai{t)^ . Inserting this into the 
Skyrme Lagrangian, one obtains the kinetic contribution to the total energy as 

T = ^a.iUijaj - aiWijhj + ^hVijbj , (31) 

where 

aj = ~i Tr TjAi^Ai , bj = i Tr rji2^2''' , (32) 
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and the inertia tensors C/y , Wij and Vij, expressed as functionals of C/o(x), are given by |17| : 



= - J Ti- i^T,T, + -[Rk,T,][Rk,Tj]j d^x, (33) 

Wi, = J ejim XI Tr [t^R^ + T;][i?fc, i?™]^ d^x , (34) 

Vij = —J eUm^jnpXlXnT'r ^RjnRp + -[Rk, Rm][Rk, Rp]^ d^ X , (35) 

where Rk = [dkUa)!!^^ is the right invariant su(2) current defined previously and 

T,= '-[n,U„]U^^ (36) 



is also an su(2) current. The potential energy, in terms of collective coordinates, is just a constant, 
the static mass of the Skyrmion. 

The quantized momenta corresponding to hi and become the body-fixed spin and isospin 
angular momenta Li and Ki satisfying the 5u(2) commutation relations 6 . The usual space-fixed 
spin and isospin angular momenta Ji and li are related to the body-fixed operators by 

.h^-D{A2)lL,, h^-D{Ai),,Kj. (37) 

We also have J'^ = L^, = K-^. Thus the operators J, L, I and K form the Lie algebra of 
S0{4)j^L ® SO{4:)j^K- Their action on Ai and A2 is given by 

[J,,A2] = ^A^n , [I^,Al] = -^T,Ai , (38) 

[L„ = -lnA2 , [K„Ai] = ^Am , (39) 

with all other commutators between momenta and coordinates zero. 

A basis for the Hilbert space of states is given by | J, J3, L3) (g) \I, I^.K^), with — J < J3, L3 < J 
and — / < /a, K3 < I. Concretely, |J, JsjLs) and |/, /aj/sTs) are Wigner functions of the Euler 
angles parametrising A2 and Ai respectively. The ground states are the states with the lowest 
values of J and / that are compatible with the FR constraints arising from the symmetries of 
the given Skyrmion. The allowed values of L3 and K3 are also constrained by the symmetries of 
the Skyrmion. In what follows, the arbitrary third components of the space and isospace angular 
momenta J3 and I3 are omitted. 

Recall that physical states satisfy [T^ 

g2^m.L|^^ ^ g27rm.K|^<, ^ (-1)^1^') , (40) 

SO even B implies that the spin and isospin, J and /, are integral, and odd B implies that they 
are half-integral. 

For general B > 1, the moments of inertia are larger for rotations than for isorotations. Since 
these appear in the denominator of the quantum Hamiltonian, the quantum states of lowest energy 
are those with minimal isospin, and spin excitations of Skyrmions require less energy than isospin 
excitations. In particular, for even B the ground state has zero isospin, but (because of the FR 
constraints) not necessarily zero spin. These observations match what is seen experimentally for a 
large range of nuclei up to baryon number about 40. 
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6 Quantization of the B = 6 Skyrmion 



The minimal energy B = 6 Skyrmion has D^d symmetry, and is well approximated by the rational 
map ansatz. The optimised rational map, in a convenient orientation, is ^1] 

Ri^)-^0^y -0.16.. (41) 

The symmetry can be seen by considering the expression for the baryon density in 114|l . The 
baryon density vanishes where the Wronskian W of the map R = p/q, given by 

W ^qp' - pq' = -2z(az* + (Sa^ - 1)2"^ + a) , (42) 

vanishes. In addition to the nine zeros of W, the baryon density also vanishes at z = c» due to the 
factor of in the denominator of the rational map. The polyhedron associated with the B = 6 
Skyrmion consists of two halves, each formed from a square with pentagons attached to all four 
sides. To join these two halves, the two squares must be parallel, with the bottom one rotated by 
7r/4 relative to the top one. The ten face centres of the polyhedron are the zeros of the baryon 
density, and the edges and vertices of the polyhedron is where the baryon density is concentrated. 

The I?4 subgroup is generated by two elements, a tt rotation about the xi-axis, and a tt 
rotation about the xi = X2 axis. The product of these is a 7r/2 rotation about the xa-axis. 
The corresponding symmetries of the rational map are R{z) — l/i?(l/z) and R{z) — —l/R{i/z). 
Therefore, the solution is invariant under the following two symmetries: 

• a TT rotation about the a;i-axis combined with a tt isorotation about the 1-axis, and 

• a TT rotation about the xi = X2 axis combined with a tt isorotation about the 2-axis. 

For the first symmetry, we have rii = n2 = (1, 0, 0)-^, and 61 — 62 — tt. As required, i?(z_n2) = 
i?ni = 1- We compute N — lb using the formula from the previous section, and deduce that 
Xfr — —1- For the second symmetry, we have rii — (0, 1,0)^, 1x2 ~ '^{^i l:0)"^i and Oi = 62 = tt. 
As required we have R{z-n2) — Rm = i- Again N = 15, so xfr — —1- Therefore the FR 
constraints reduce to 

e*''^ie"-^i|*) = -I*) , (43) 

and 

gi7r(Li+L2)/%/2gi7rK2|^^ ^ _|^^ ^ (•44') 

We note that these constraints are equivalent to those obtained in 15^ using an alternative 
method by which the contractibilities of the closed loops corresponding to the symmetry group 
elements were determined by continuous deformation of the minimal energy solution into three 
well-separated B — 2 tori. 

A basis for the Hilbert space of states is given by |J, L3) (81 lIjK^). We require the identities 
HHI: 

e-^H-/,i3) = ^|J,4)i?4^3(0,^,7r) = (-l)-V,-i3), (45) 
e"^M^,^3) = ^\I,K',)D'^,^^iO,7:,n) = {-iy\I,-K,), (46) 
e-^=|/,X3> = ^|/,if^)i?i,^3(0,^,0) = (-l)^+^^|/,-if3>, (47) 
"'^^+^^^/^l-/,i3) = ^|J,4)^4^3(0,^,7r/2) = (-l)'^+^^^'|J,-i3), (48) 

L' 



e 
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where Dj^, , for given J, is the matrix of Wigner functions representing the symmetry operation. 

Seeking simultaneous solutions of the above FR constraints, we obtain the ground state as 
|1, 0) (g) |0, 0). The first excited state with / = is |3, 0) (g) |0, 0), and the lowest state with / = 1 is 
|0,0)(^|1,0) US]. 



7 Mean Charge Radius of the Quantized B = 6 Skyrmion 

The state |1,0) |0,0) of the B = 6 Skyrmion has the quantum numbers of the ground state of 
the Lithium-6 nucleus. Let us calculate its mean charge radius (r^)2, defined as the square root 
of 

„2\ _ ^ /,Trl / „2 7 /„ ,\ j3 



r- = 



^ J r^Joi^,t)d''x\-^), (49) 



where Jo is the electric charge density operator and represents the state with any value of J3. 
Only the isoscalar part of jo(x, t) contributes to the matrix elements, and so the integrals are pure 
c numbers. The normalization factor is 3, the total electric charge of Lithium-6. Therefore 

{r^)^\jr^B^{^)d?x, (50) 

using |(SJ). 

Using the expression (|14|l for the baryon density in terms of the rational map ansatz, we obtain 

= I /"(-/'jr'siiiV*. (52) 

The profile function / can be numerically determined using a shooting method^, for which the 
inputs are B, m and I. Then the radial integral in (j52() can be evaluated, and the mean charge 
radius in Skyrme units determined as the square root of this. B — 6 and T — 50.76 for the rational 
map (|41() . but we defer discussion of m and the actual value of (r^)2 until section |5| 



8 Quadrupole Moment of the Quantized B = 6 Skyrmion 

An additional static electromagnetic property that is known experimentally for the Lithium-6 
nucleus to quite high precision is its quadrupole moment. The classical quadrupole tensor of a 
Skyrmion is defined as 

Qab = J {iXaXb ~ r'^Sab)Jo{^, t) X = ^ J i^^aXb - r^(5afc)Bo(x) X , (53) 

using ((SJ. It is traceless by definition. The symmetry of the B = 6 Skyrmion implies that the 
tensor Qab is diagonal, and further that Qn = Q22 = ^Qss/S. 

From H53(l we find 

Q33 = IJ r^{icosH-l)Bo{^)d^x, (54) 

^Thanks to Bernard M. A. G. Piette, University of Durham, UK, for providing the CH — h code 
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and therefore, using the rational map ansatz, 



Q33 = i^/(3C0S^.-1) (1 



\R\' 



dR 



dz 



sin 0d0d(l) - I (-/') sin^ / dr . (55) 

TT 



For the rational map given previously for the B = 6 Skyrmion, 



1 + 



dR 



dz 



A\zmi + \z\'y 



(az« + (3o2 - l)z4 + a) {dz^ + (Sa^ - l)^^ + a) 
(|a|2|zP + |z|8 + |z|4(l + az* + dz^) + dz^ + az* + \a\^f 



(56) 



and substituting z — tan(^?/2)e"^, writing a = ia, where a G M, and setting p — 1 + 3a^, this 
becomes 



1 + 



dR 



dz 



= 4 tan"' 



1 + tan^ 



a^tan^^ f + 2Q;^sin40tani2 f + (/^^ + 2a'^ cos8(j)) tan*^ f - 2a/3sin40tan 



4« 



(57) 

[ (a^tan^^f + (l-2asin40)tan^f + (l + 2asin4(/i)tan''f +a2) J 

A numerical integration technique can then be used to determine the angular integral in H55|l . for 

2 



a = 0.16. The result is 



0.395 ■ 



(-/')r2 sinVdr-, 



(58) 



where the radial integral is the squared charge radius 1)52(1 . The accuracy of the numerical integra- 
tion can be checked by considering the integral 



1 



l + |i?P 



dR 



dz 



sin d0 ( 



(59) 



The same procedure yields a result 75.40, which is equal to 247r. 



Note that Qsa is positive, so the classical Skyrmion is prolate. It is helpful to consider upper 
and lower bounds on its value, for a given B and (r^). We observe that is maximal for baryon 
densities restricted to the xa-axis and minimal for baryon densities restricted to the (xi, a;2)-plane, 
hence 

-|(r^)< 033 <S(r^). (60) 



For the B = 6 case, 

-3<Q33/(r-') <6. 
Within this range, the actual value Q33/ (r^) — 0.395 is rather close to zero. 

In the quantum state the Skyrmion occurs in all possible orientations. Substituting 

J7(x) = AiC/o(i?(.42)x)4, 

we obtain 



Q,, - D{A2ZQabD{A2) 



(61) 



(62) 



(63) 



for the rotated classical Skyrmion. We observe that Qij is independent of Ai. The symmetry 
relations reduce this quadrupole moment operator to the expression 



Q^, = ^Q33 {3D{A-2)lD{A2)3, ~ 6,,) , 



(64) 
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and thus 

Q33 = ^Q33 (3 COS^ 0-l)= Q33i^oo(</', ^) : (65) 

where 0, 0, ip now denote the Euler angles parametrising A2, and I?gQ(0, 9, ip) is a Wigncr function. 

The wavefunction 5* is the product of Wigner functions D^^^cf), 9, ip)®Dl^{a, /3, 7) on 50(3)'^ x 
SO['iY . For the ground state of the 5 = 6 Skyrmion, with J = 1 and 1 = 0, and J3 = m, 



(66) 



The quadrupole moment of the quantized _B = 6 Skyrmion is defined as the expectation value of 
(333 in the state with J3 = 1 (by convention the quadrupole moment is measured in the top spin 
state), i.e. 

Q = (^'j3=i|Q33|*J3=i) - ^Q33 j Dlo\Dl,\^ Sin 9 del) d9d^, (67) 

and using the well-known expression in terms of Wigner 3j-symbols for the integral of the product 
of three Wigner functions, 

J f j" \ ( J f j" 



Di,Di,Di^sm9d^d9d^P = 8n'^'^ e j (, & d / j' ^^^^ 

Q33 J Dl„Dl,Dl_, sin 9 d9d(l>d^, as 0^,^ ^ {-ir'-^'DL 



we find that 

■ 2 1 1 \ / 2 1 1 

j I 1 -1 



= -sg 



33 



= -^^33. (69) 



This is the standard result for an axisymmetric system in a J = 1 state. Note that a prolate 
classical shape {Q33 positive) gives an oblate quadrupole moment {Q negative), and vice versa, 
because the prolate object can be regarded as spinning about its short axis aligned along the third 
axis in space when J3 = 1. From ()61|l we deduce the inequalities 

-1.2 < Q/(r2) < 0.6. (70) 



We note that to obtain physically meaningful results for both the mean charge radius and the 
quadrupole moment, it is necessary to convert back to physical units, by reintroducing the Skyrme 
model parameters, and using the conversion factor he — 197.3 MeVfm. We do this in the next 
section. 



9 Reparametrising the Skyrme model 



The parameters e and F^^ can be fixed in a number of ways. It has been common practice to use 
the set of parameters given in specifically 

e = 4.84, = 108 MeV and = 138 MeV (which implies m = 0.528) (71) 



11 



for the Skyrme model with the physical pion mass taken into account. In the values of e and 
were tuned to reproduce the masses of the proton and the delta resonance. This parameter 
set was adjusted to optimise the predictions of the model in the B — 1 sector at the expense of 
the B = sector, which requires = 186 MeV. It is therefore not the optimal parameter set 
with which to describe the higher baryon number sectors of the Skyrme model. Nevertheless, in 
[H] this parameter set was used in a quantization of the collective coordinates of the toroidal B = 2 
Skyrmion. The resulting quantum state is too small compared to the size of the deuteron, and 
it is too tightly bound. However, this discrepancy was dealt with by the quantization of selected 
vibrational modes of the B = 2 Skyrmion in . We therefore do not propose to change e and 
for B = 1 or i? = 2. However, for higher B, it becomes technically very difficult to perform a 
quantization of the possible vibrational modes. Instead, we consider rescaling the parameters in 
order to fit the Skyrme model to experimental nuclear data, whilst suppressing the quantization 
of vibrational modes. This is a type of renormalization. 

We have first determined (r^)5 and Q for the B = 6 Skyrmion using the set of parameters 
given in [T], and setting a = 0.16i (the value of a that minimises 2, as required by the rational 
map ansatz). After evaluating the radial integral l|52(l . we find 

(r^) 5 = 1.48 fm, Q = -0.173 (fm)2 . (72) 

In |1U| the experimental root mean square radius of the charge density distribution of the 
Lithium-6 nucleus is given as (r^) 2 = 2.55 fm, substantially larger than 1.48 fm, and its quadrupole 
moment is given as —0.82 x 10"'^ barns — —0.082 (fm)^ in (HI. We can change the length scale 
2/ei^7r, to obtain a mean charge radius in close agreement with experiment. This results in the 
dimensionless mass m — 2m7r/ei^7r, with 771^^ — 138 MeV as before, being changed, and so the 
profile function and hence the radial integral in the expressions for the mean charge radius and 
quadrupole moment is also modified. By iteration, it is found that by setting m — 1.125, and 
hence eF^r = 245.34 MeV, we obtain (r^)2 = 2.55 fm, which fits the data. 

With this choice for eF^, the theoretical prediction of the quadrupole moment becomes Q = 
— 0.512 (fm)^. Although the sign is right, this is unfortunately about six times the experimental 
value. One must bear in mind that the rational map ansatz does not generally determine exact 
solutions of the Skyrme field equation, but rather approximations to solutions. It is useful to allow 
a slight modification of the rational map l|41() . while preserving its symmetry and prolateness. If 
we take a = 0.1933i, rather than 0.16i, and use the new value for eF^^, we obtain a quadrupole 
moment in very close agreement with experiment: Q = —0.082 (fm)^. 

We note that this modification of the rational map leads to the value of the integral T being 
slightly modified (from 50.76 to 51.49, i.e. not quite the minimum as required by the rational 
map ansatz) , which subsequently leads to the profile function and hence the radial integral in both 
the expressions for the mean charge radius and the quadrupole moment being modified, but only 
slightly. In summary, using eF^^ = 245.34 MeV, ttItt = 138 MeV (which implies m = 1.125), and 
a = 0.1933j, we find 

(r^) 5 = 2.55 fm , Q = -0.082 (fm)^ , (73) 

and hence Q/{r'^) = —0.0127. By virtue of 1)70(1 . we consider the quadrupole moment of the 
Lithium-6 nucleus to be only slightly oblate. The change we needed to make to the parameter a 
should not be regarded as a large change. 

It remains to consider the energy scale FT^/'ie of the Skyrme model. Assuming negligible spin 
kinetic energy (it is expected to be in the region of IMeV), we may identify the static Skyrmion 
energy 1(18(1 with the mass of the Lithium-6 nucleus, Af, given by 

M ^ QmN - EunMng = 5601 MeV, (74) 
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where ttin = 939 MeV is the average mass of a nucleon and Etinding = 32 MeV is the total binding 
energy of the nucleus |23j . The static Skyrmion energy is determined numerically, and we find for 
m — 1.125, after converting to physical units, that 

= . 972 . (75) 

4e ^ ' 

Equating this to the mass of the nucleus implies that 

i^^/4e = 5.76 MeV, (76) 
and combined with eF^r = 245.34 MeV we obtain, finally 

e = 3.263, = 75.20 MeV and = 138 MeV (which implies m = 1.125) . (77) 



We conclude that with this choice of parameters, and taking a — 0.1933i, the mass and static 
electromagnetic properties of the quantized B = 6 Skyrmion are in close agreement with those 
experimentally determined for the Lithium-6 nucleus. As mentioned previously, this value of F^^ is 
considered to be a renormalized pion decay constant, and the fact that it is so much less than the 
experimental value is therefore not necessarily conflicting. 



10 Reconsidering the B = 2 and B = 4 Skyrmions 
10.1 The quantized toroidal B = 2 Skyrmion 

The minimal energy solution in the B — 2 sector is axially symmetric, and has the shape of a 
torus. In 1^, Braaten and Carson performed the collective coordinate quantization of the solution 
(see also [El)- Using the same notation as that used above for the B — 6 case, it was found that 
the axial symmetry of the solution requires 

(2if3 + i3)|*) =0. (78) 

In addition to this constraint, the FR constraint associated with the discrete rotations by tt in the 
D2 symmetry group is 

e^'^^^e'^^^l^') = -I*). (79) 

The ground state is the J = 1, / = state |1,0) (g) |0, 0), which has the quantum numbers of the 
deuteron. The first excited state |0, 0) ® |1, 0) may be identified with the isovector ^5*0 state of the 
two-nucleon system. 

Further properties of the quantized Skyrmion were also determined, taking the Skyrme model 
parameters to be as in (|71|l . and compared to experimental data on the deuteron. An expression 
equivalent to H5(J|I is obtained for the mean charge radius of the _B = 2 Skyrmion. The expression 
for the quantum quadrupole moment in this case is also given hy Q — — -gQas, where Q33 is negative 
here, because the classical Skyrmion is oblate. The quantity Q33 was evaluated in using the 
numerical, exact classical solution and its corresponding baryon number density. The predictions 
were: 

(r2)5 = 0.92 fm, (80) 
Q = 0.082 (fm)2, (81) 
M = 1720 MeV, (82) 
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compared to experimental data of (r^) 5 = 2.095 fm[TT|, Q = 0.2859 (fm)^ ^ and M = 1876 MeV. 
The mass of the nucleus (I82f) is given by the sum of the classical mass, equal to 1659 MeV, and 
a spin correction equal to l/Vn — 61 MeV. It was suggested in 6 that these discrepancies could 
be an artifact of the Skyrme model parameters, and that the disagreement with experiment could 
be assuaged by an adjustment of e and -FV such that C-FV is approximately halved and -FV/e is 
unchanged. The new parameter set which we propose is close to this, but we also take into account 
the corresponding change of m when fitting to experimental data. 

10.2 Attractive channel Skyrmions and the deuteron 

Because the deuteron appears too small and too tightly bound according to 6 , it was suggested 
in |19| that the deuteron is better described as a quantum state on a ten-dimensional manifold of 
B = 2 Skyrme fields. This manifold includes the toroidal configurations of minimal energy (i.e. 
the eight-dimensional manifold considered in |^), and deformations of these into configurations 
which are approximately the product of two separated B — 1 Skyrmions in the most attractive 
relative orientation. The ten coordinates of this manifold are the separation parameter, overall 
translations, rotations and isorotations. A unique bound quantum state is found with the quantum 
numbers of the deuteron. The values of e and F^^ were taken to be 4.84 and 108 MeV respectively, 
i.e. those of ^ and 0, but for technical reasons the pion mass was taken to be zero. In this model, 
the predictions are 

(r^)^ = 2.18 fm, (83) 
Q = 0.83(fm)^ (84) 
M = 1872 MeV. (85) 

The calculated deuteron binding energy of 6 MeV is very close to the experimentally determined 
binding energy of the deuteron, but the total mass could not be straightforwardly determined. 
The prediction of the deuteron mass H85|) needs to be interpreted as that of the sum of the proton 
and neutron masses minus this binding energy. The prediction for the mean charge radius is in 
very close agreement with experiment. However, the deuteron's quadrupole moment has become 
much larger than the experimental value. This is explained in (19) as follows: The tensor force is 
responsible for the existence of a d-wave contribution to the deuteron wave function. In the absence 
of a d-wave, Q would be zero. The size of Q is therefore a measure of the d-wave probability, which 
in turn indicates the strength of the tensor force. The truncation to the space of attractive channel 
fields systematically overestimates the strength of the tensor force. Therefore it is not surprising 
that the theoretical prediction for Q is rather large. 

10.3 Is reparametrisation helpful in the B = 2 sector? 

Let us reconsider the collective coordinate quantization of 0, ignoring vibrational modes, but 
using our new parameter set (|77|) . To avoid recalculating the exact B — 2 Skyrmion solution for 
various values of the pion mass parameter m, we use the rational map ansatz. There is a unique 
map here, namely R{z) ~ |14| . This gives, for the quantized B — 2 Skyrmion, 

^ Ji-Dr' sinVdr, 

(86) 

f')r^sm^fdr. (87) 



1 1 

2 4^ 



1 + 



dR 

dz 



sinddOdcj) - j {-f)r^ sin^ f dr 



and 



g^4i^/(3cos- 



dR 



dz 



sin eded(j} - / (- 

TT 
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It is straightforward to calculate the profile function / in this case, and to perform the numerical 
integrations. Converting back to physical units, and using first the traditional values of e, F^^ and 
m given by H71|) . we obtain 

(r2)5 = 0.940 fm, (88) 
Q = 0.102 (fm)^ (89) 
M = 1757 MeV. (90) 

which are reasonably close to the values (|80|I ~ H82|I predicted in We have again assumed that 
the mass of the nucleus H90|l can be equated to the sum of the static Skyrmion energy (1696 MeV) 
and the spin correction which we take to be the same as in jSj (61 MeV). Comparing this value of 
the static Skyrmion energy to the classical mass given in [H], we see that the rational map ansatz 
gives a very good approximation. 

Using the new parameters (|77|l . and their associated length and energy scales, we find after 
recalculating the radial integral (which we recall is a function of m) 

(,-2^5 ^ 1.644 fm, (91) 
Q = 0.311 (fm)^ (92) 
M = 1969 MeV. (93) 

The deuteron mass (|93|l is again assumed to be equal to the static Skyrmion energy (1950 MeV) plus 
an estimated spin correction equal to 61 MeV times the ratio of the squares of the predicted mean 
charge radii (|80|l and (|91l) , which allows for the increase in the moment of inertia. Comparing these 
with the experimental values, we conclude that the new parameters give quite a good description 
of the deuteron, better than those in [B], and provide an alternative to taking account of the soft 
vibrations. 



10.4 The B = A Skyrmion and the a-particle 

The B — A Skyrmion is cubically symmetric, and was first quantized by Walhout |26|, using 
the exact classical solution and the traditional parameter set. The rational map ansatz for this 
Skyrmion uses the unique map of degree 4 with cubic symmetry |14j . 
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^ + 2V3^z , . _ 
^ ' z4_2V3iz2 + l ^ ^ 

The FR constraints corresponding to the cubic symmetry were determined in |15j as 

The ground state, in the basis introduced previously, is |0,0) (E) |0, 0), which trivially satisfies the 
constraints. This state may be identified with the a-particle. 

If vibrational modes are ignored, then the traditional parameters give a much too small and 
tightly bound a-particle. Walhout's rather complicated analysis takes into account a number of 
the vibrational modes of the B ~ A Skyrmion, obtaining {r'^)2 = 1.58 fm and M = 3677 MeV, 
compared with the experimental values {r'^)^ = 1.71 fm JOI and M — 3727 MeV. 

We wish to ignore the vibrational modes, but in compensation, use our new parameter set H77f) . 
As the spin of the B — A Skyrmion is zero, there is no quadrupole moment to determine, just the 
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mean charge radius. The mass of the quantized Skyrmion is equal to the classical mass, as there 
is no spin energy contribution. Numerically computing the profile function (with inputs X — 20.65 
and m — 1.125), and numerically integrating the corresponding densities, we find 



which agrees reasonably well with experiment. 



11 Conclusion 

We conclude that in order to accurately predict certain properties of nuclei (with baryon numbers 
greater than two) using the Skyrme model, it is necessary to rescale the Skyrme model parameters, 
if the quantization of vibrational modes is not going to be considered in the calculation. We have 
in this paper explicitly determined the values of these parameters in the B = 6 case, and discussed 
how well they work in the B = 2 and B — 4 cases. It would be interesting to see how well they 
work for Skyrmions with other baryon numbers. 

Our suggestion that the dimensionless pion mass m should be increased from 0.528 to 1.125 
backs up the suggestion made in in which it was found that as the pion mass is increased 
its effect becomes important for Skyrmions of all baryon numbers. In particular, the hollow, 
polyhedral shell-like solutions which exist for a large range of baryon numbers up to 20 or 30, and 
are stable for zero pion mass, become unstable for baryon number eight and above when m is of 
order unity. Instead, the stable solutions become more solid structures, some of which are related 
to chunks of the Skyrme crystal. This improves the qualitative fit of the Skyrme model to real 
nuclei. Also, in low energy Skyrmion solutions composed of charge four sub-units were found 
for baryon numbers a multiple of four, with m = 1. This is the Skyrmion analogue of the a-particle 
model of nuclei. 

The standard values of the Skyrme parameters originate from in which attention was 
restricted to the single nucleon. In recent years, the substantial increase in results available for 
Skyrmions over a range of baryon numbers enables us to fit the Skyrme parameters to experimental 
data for larger nuclei, with a view to modelling nuclei of all baryon numbers. The parameters we 
propose here should be regarded as provisional, since they rely on just a few basic properties of 
Lithium-6, and depend on the rational map approximation. More work remains to be done within 
the Skyrme model to see if one set of Skyrme parameters can fit experimental data over a range of 
baryon numbers. Both static properties of nuclei, like their masses and electric charge distribution, 
and also the excitation energies of higher spin states, need to be considered further. 
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(r2)2 ^ 2.13 fm, 
M = 3679 MeV , 



(97) 
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